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We show model wavefunctions for neutral collective modes in fractional quantum Hall (FQH) 
states have simple analytic forms obtained from judicially reducing the powers of selected pairs in the 
ground state Jastrow factor. This scheme of "pair excitations" works for the magneto-roton modes 
of single-component Abelian and non-Abeliean FQH states, as well as neutral fermion mode for the 
Moore-Read (MR) state. The analytic wavefunctions enable computations in the thermodynamic 
limit previous inaccessible to numerics, and the long wavelength limit of the neutral energy gap of 
the magneto-roton modes can be interpreted as fusion of charges in two-dimensional plasma picture, 
extending the plasma analogy to neutral excitations. A lattice diagrammatic method of representing 
these many-body wavefunctions and FQH elementary excitations is also presented. 



The fractional quantum Hall effect (FQHE)[T] is one 
of the prime examples where strong interaction between 
electrons dictates the dynamics. A very fruitful approach 
for such a strongly correlated system is to look for model 
wavefunctions and model Hamiltonians, where physical 
interpretations are more transparent, that are adiabat- 
ically connected to experimentally accessible systems. 
For the ground states and charged excitations of FQHE 
at several filling factors, wavefunctions can be written 
down analytically [14- I6j. We can thus identify the prop- 
erties of FQHE in the thermodynamic limit, and rein- 
terpret the wavefunctions in analogy to two-dimensional 
plasmas [14], or as conformal blocks of conformal field 
theory (CFT) [HI [26] . In general for these approaches, 
incompressibility of FQHE is always assumed, and the 
dynamics is not explicitly discussed. 

Incompressibility of FQHE is defined by the neutral 
collective excitations. Such excitations are important for 
understanding which topological phases of FQHE can 
be stabilized, and also very relevant to recent develop- 
ment of fractional Chern insulator, where the issue of in- 
compressibility and its mapping to FQHE [3 -9] are areas 
of active research. The first formal treatment of neu- 
tral excitations came from single mode approximation 
(SMA) for the magneto-roton[2] mode, where good model 
wavefunctions of density wave excitations can be con- 
structed numerically up to the momentum of roton min- 
imum, beyond which SMA is no longer valid. Following 
experimental studies of the collective modes by several 
groups [HUH 3] , more recently the model wavefunctions 
are constructed numerically both for the magneto-roton 
modes in the Read-Rezayi (RR) series, and the neutral 
fermion mode[T7] in the Moore-Read (MR) state charac- 
teristic of its non-Abelian nature. These wavefunctions 
are good for the entire range of momentum. However in 
practice the long wavelength limit is not accessible due 
to limitation by the system size. It is now understood 
that the collective modes can be thought as excitons of 
composite fermions[20 22], and in the long wavelength 
limit it is a "spin 2" quadrupole excitation; as momen- 



tum increases it relaxes into a dipole excitation beyond 
the roton minimum[18] [23]. Interestingly, even though 
the underlying phenomenological pictures of the collec- 
tive modes can be different, it is found numerically that 
exactly the same set of model wavefunctions are found 
with different approaches, with very rich algebraic struc- 
tures. This suggests a natural way of representing these 
model wavefunctions for the neutral excitations as well, 
with no need of explicit variational parameters. 

In this Letter, we present analytic wavefunctions that 
are identical to those numerically generated in [2"Ull23| . 
and calculate the energy gap of the quadrupole exci- 
tations in the thermodynamic limit. We start by pre- 
senting the wavefunctions of the collective modes for 
fermionic Laughlin state at filling factor v — 1/m in 
the lowest Landau level (LLL), where m is odd. On 
the sphere the ground state is the Laughlin wavefunc- 
tion in total angular momentum L = sector, or the 
fermionic Jack polynomial Jj^oiooi |19]. By stripping 
away the single particle normalization factor the holo- 
morphic part of the wavefunction is the same on all 0- 
genus manifold. Thus we label the state with its to- 
tal angular momentum on the sphere. Taking z = 
7^7 ( x + i y)' z * = ^W^( x ~ i y)' where 1 b = y/h/eB is 
the magnetic length, the unique ground state is given by 
ipl = Yli<j( z i — Zj) m e~2 ZiZi , with model Hamiltonian 
made of Haldane pseudopotentials [2"4] V — J2i<jViji 
with 

/i2 j2 rn—1 
^^i n ( 9 2 /|) e -^ 2 ^e^^-^) (1) 

71 = 

where L n {x) is the n th Laguerre polynomial and Ri is 
the guiding center coordinate of the i th particle. Physi- 
cally, Vij is the short range interaction that projects into 
the two-body Hilbert space with their relative angular 
momentum smaller than m. The family of collective 
modes at different angular momentum sector (we omit 
the exponential part of the wavefunction, which is ir- 
relavent in LLL) is as follows: 



2 



A[{z 1 -z 2 r' 2 \\(z l -z j r\ l = 2 

i<3 

l 

A[(zi - z 2 r~ 2 ( Zl - zs)™- 1 Y[(zi - z 3 ) m ] L = 3 

i<j 

l 

A[(zi - Z2) m - 2 ( Z1 - z s ) m - l (zi - zi)™- 1 H(zi - Zj ) m ] L = 4 




FIG. 1. For N e particles, the lattice can be viewed as an 
N-gon, with three bonds connecting every pair of vertices 

Here A indicates antisymmetrization over all particle 
indices, and Yii<j means products of only pairs {ij} that 
do not appear in the prefactors to the left of it. Thus the 
L = 2 state, which is the quadrupole excitation in the 
thermodynamic limit [22], is obtained from the ground 
state by reducing the power of one pair of particles (which 
we choose arbitrarily as particle 1 and 2) by two, followed 
by antisymmetrizing over all particles. This scheme nat- 
urally forbids an L = 1 state by pair excitation, since 
if we reduce the power of one pair of particles by one, 
antisymmetrization kills the state. 

The L = 3 state is generated by pairing particle 1 with 
another particle and reducing their pair power by one. It 
is now clear how the modes in other momentum sectors 
are generated. Naturally for a total of N e particles, the 
family of collective modes ends at L = N e , agreeing with 
the scheme in [23 . Indeed all wavcfunctions here sat- 
isfy the highest weight condition, and the states relax to 
the ground state far away from the excited pairs; these 
are exactly the conditions we used to numerically gener- 
ate the unique model wavefunction in each momentum 
sector. 

An intuitive way to visualize the family of collective 
modes is to map the particles onto a lattice, where each 
lattice site represents a particle. Since for FQHE we have 
a quantum fluid instead of a solid, every two lattice sites 
interact with each other. The number of bonds between 
each pair of lattice sites equal to the power of the pair of 
particles in the wavefunction. As an example we consider 
the simpliest Laughlin state at v = 1/3, so for the ground 
state every two lattice sites are connected by three bonds, 
as shown in Fig. [TJ 

The collective modes are obtained by breaking the 
bonds between lattice sites, as shown in Fig. [2j We can 
view the entire family of collective modes as elementary 
excitations centered around a single red lattice site. Note 
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FIG. 2. (Color Online) Collective modes from L — 2 to L — 5, 
where the change of bonds are highlighted with red color. 



the lattice pattern uniquely defines the many-body wave- 
function, and different types of "elementary excitations" 
can be identified with different patterns of bond-breaking 
around a single lattice site. 

The same scheme applies to MR state. It is instruc- 
tive to first see how the MR ground state is obtained. 
The Laughlin wavefunction at half filling is given by 

^ioioioi— ( z i) ~ Yli<j( z i~ z j) 2 ■ F° r fermions this is not a 
valid state; instead the ground state was constructed by 
a pairing mechanism|15 , which is also a Jack polynomial 
^iiooiiooii- - • The pairing reduces the power of each pair 
of particles by one. Explicitly we have for 2n particles 



Y[( z i ~ z j) 2 ~* - Z 2)(Z3 - 24) ■ ■ ■ (Z2n-1 - Z2„) Y[( z i _ z i) 2 \ 

=p f (_i-)n(*-*i) a (3) 

where the last line of Eq. ([3| is the familiar Pfaffian for 
the MR ground state. This suggests lattice representa- 
tion of MR state and its magneto-roton mode with the 
same scheme, as shown in Fig. [3] 

We can also use the same scheme to generate the neu- 
tral fcrmion mode for the MR states with odd number of 
particles. In this case, starting from the Bosonic Laugh- 
lin wavefunction at half filling, every two particles form 
a pair except for one particle. Naturally the "ground 
state" of the neutral fermion mode is given by 
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FIG. 3. (Color Online) The lattice configuration of the ground 
state L = and the first collective mode L = 2. Consecutive 
collective modes can be obtained by breaking one of the double 
bonds connecting the red lattice site to some other site. 




For even number of electrons we have N e = 2n and 
Eq.(|6]) is the MR ground state. The magneto-roton 
modes are given by 

Vw fe+2 = Y[(zi - ZjfAlV^VlaVis ■ ■ ■ Pi.a+2*] (7) 

For odd number of electrons we have N e = 2n + 1 and 
Eq.(|6]) is the MR quasihole state of Eq.Q. The neutral 
fermion modes are given by 

_3 N " 

Vm7 5+k = Y\(z t - z 3 ) 2 A[P^P 2 NcA r N<! , 3 ■ ■ ■Vn^i+u] (8) 

i<j 

The analytic wavefunction is useful in calculating the 
magneto-roton mode energy gap in the long wavelength 
limit. For the Laughlin state, the energy gap is given by 



FIG. 4. (Color Online) The lattice configuration of the zero 
mode quasihole state L = ~ (N e — 1) and the first two neutral 
fermion modes at L — | and L = §. Consecutive collective 
modes can be obtained by breaking one of the double bonds 
connecting the red lattice site to some other site. 



Y\i Z i ~ Z if ~^ - Z 2)(Z3 - Z&) ■ ■ ■ (Z2n-1 - Z 2n ) 
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For model three-body Hamiltonian, this is a zero- 
energy abelian quasihole state Juoonooii---ooiiooi m ^ ne 
angular momentum sector L = \{N e — 1). Its lattice pre- 
sentation and that of the neutral fermion mode is given 
in Fig. [4] 

To write down the analytic wavefunctions in a more 
formal way, we define Vij = z ._ z . ■ Notice the Pfafnan 

for 2n particles can be written as Pf ^ - \ z / j ~ .A['P < ' 2 "* 1 ], 

where V^ 2 ^ = V\iV^ ■ ■ ■ V2n-i,2n- The magneto-roton 
mode for the Laughlin state is given by 



We already know from [23] that in L = 2 and L = 3 
sector, SMA is exact for the magneto-roton mode model 
wavefunction. Defining the guiding center ladder opera- 
tors as b\ — Zi,bi — d Zi , we have 

V^ 2 = 9 f 1 tt E^) 2 ^ (10) 
2m (m — 1) * — ' 

i 

In the thermodynamic limit, the normalization con- 
stant of the above two modes are thus related to the long 
wavelength expansion of the ground state guiding center 
structure factor: 



= -^(g ab qaq b ) 2 + 0(q e ) (11) 

where s = is the guiding center spin0 |2"3] . and 
g ab is the guiding center metric [25] . We thus have 
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where S is the symmetrization over all particle indices. 
From the Bosonic Laughlin wavefunction at filling factor 
1/2 we can impose pairing to obtain 
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The numerator of Eq. ^ can be calculated with plasma 
analogy. Note in Eq.Q, before antisymmetrization the 
term only has one pair of particles with relative angular 
momentum smaller than m. We thus have 
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-mvl^VlMi) (13) 
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where Af is the normalization constant of the Laughlin 
state. We note that V12 projects out states with relative 
angular momentum (z%— Z2) m ~ 2 , which can be integrated 
over. The numerator is thus equivalent to evaluating the 
norm of the following wavefunction: 



n 



V2 



n 

Ki<j<N e 



{zt-ZiT (14) 



which can be evaluated as the free energy of two- 
dimensional one-component plasma (OCP) on a disk 
with radius R 2 — and elementary charge e = 

2\/irmkBT 1 where particle 1 interacts with the rest of 
the particles with charge 2e. We thus obtain 



l m(m — l) 2 ^2-^ 

L_ k B T 

7TS 



(15) 



Both F2 and F are free energies of OCP in the thermo- 
dynamic limit (N e — ¥ 00), where F is for N e particles, 
each with charge e with logarithmic two-body interac- 
tions together with a neutralizing background of radius 
R; for J-jjj we have the same neutralizing background 
with N e — 2 particles of charge e, and exactly one parti- 
cle with charge 2e. Thus F2 — F is the free energy cost 
of fusing two particles of charge e to create a particle of 
charge 2e, which is an 0(1) effect. 

Similar calculations can be carried out for the 
magneto-roton mode in the MR state. Analogous to 
Eq.(10l we have ^= 2 = §jV\6?Vw, an d in the long 
wavelength limit we have 
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where s mr = — 2 is the guiding center spin for the MR 
state, and Fu is the standard two-component plasma 
free energy for the MR ground state [55]. The charge 
for interaction between the two components is given by 
Qi = ±v / 3fcj3T, while the charge for interaction between 
one component and the neutralizing background is given 
by Q2 = 2yHtBT. F3 — Fu is the free energy cost of 
fusing three particles for each component to create one 
particle with charge 3Q2 but with the same ±Qi- 

The evaluation of the long wavelength gap of the neu- 
tral fermion mode is less transparent. The difficulty 

L=- 

lies with evaluating the normalization constant of ip mr 2 . 
There is no SMA for the neutral fermion mode, and it is 
not known if in the thermodynamic limit the gap should 

be inversely proportional to the guiding center spin. On 

l=- l=s 
the other hand (^ mr 2 1 V3bdy |Vw 2 ) can be mapped to 

2-component plasma as well, and we obtain 



Here Fu is the free energy of the 2-component plasma 
similar to that of Fu with only one difference: there is 
exactly one more particle carrying charge Q2 that inter- 
acts with the neutralizing background, and its Q\ charge 
is zero. This is how an unpaired fermion in the MR 
state is interpreted in the plasma analogy. Furthermore, 
F3 — Fu is the energy cost of fusing the unpaired fermion 
with one pair of two other fermions, creating a particle 
with charge Q2 = Q^JksT but again with zero Q\. 

In conclusion, analytic wavefunctions for both the 
magneto-roton modes and the neutral fermion modes are 
presented. The energy gap of the quadrupole excitation 
in the thermodynamic limit can be related to the free 
energy cost of the fusion of charges in the plasma en- 
ergy, and is inversely proportional to the guiding center 
spin which characterizes its topological order. This is 
the first time that the plasma analogy is extended to 
neutral excitations of FQHE, and the analogy not only 
applies to the wavefunctions, but also to the dynamics as 
well. The lattice diagrams we presented uniquely defines 
the many-body wavefunctions; one would also conjecture 
the diagrams are useful in determining the many-body 
wavefunctions of multi-roton excitations. Since the col- 
lective mode in the long wavelength limit is buried in 
the multi-roton continuum, it is important to calculate 
the decay rate of the collective mode. Numerical calcula- 
tion has been performed to show that even in the contin- 
uum the decay rate of the collective mode is very small. 
This opens up the possibility of experimental detection 
of these modes. A more detailed analysis of the decay 
rate of collective modes will be presented elsewhere [27]. 
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